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MIXED  CONDUCTION  IN  SEMI-INSULATING  GALLIUM  ARSENIDE 


INTRODUCTION 

With  the  ever-increasing  use  of  gallium  arsenide  in  a  wide  variety  of  integrated 
circuits  and  millimeter-wave/microwave  devices,  a  demand  has  arisen  for  reproducible 
semi-insulating  substrate  material  for  device  fabrication.  Synthesis  of  such  materi¬ 
al  has  been  hampered  by  an  incomplete  and  sometimes  erroneous  knowledge  of  its  trans¬ 
port  properties,  and  the  impurity  states  by  which  they  are  affected.  Especially  trou¬ 
blesome  has  been  the  difficulty  in  separating  the  contributions  of  electrons  and 
holes  in  the  range  of  Fermi  levels  where  neither  carrier  dominates  conduction,  that 
is,  the  mixed  conduction  region  (MCR).  An  early  attempt  at  a  solution  met  with  some 
success,^ but  was  found  inadequate  for  many  samples  since  it  neglected  single  car¬ 
rier  magneto-resistance  and  magneto-Hall  effects.  A  new  version  developed  a  family 
of  curves  relating  electron  mobility  to  Hall  mobility  through  resistivity.3  That 
analysis  and  those  employed  by  other  workers  rest  on  assumptions  or  semi -empirical 
deductions  about  the  electron/hole  mobility  ratio  in  conjunction  with  either  the  in¬ 
trinsic  carrier  concentration  or  the  actual  electron  mobility.  Divergent  views  of 
criteria  establishing  the  MCR  result,  and  these  are  masked  in  part  by  the  different 
temperatures  assumed  in  the  various  analyses.  However,  the  adopted  mobility  ratios 
used  in  these  analyses  are  similar, and  can  be  used  with  established  values  of  the 
energy  gap  and  the  effective  densities  of  states  at  the  band  edges  to  form  a  more 
general  and  comprehensive  analysis  applicable  at  all  reasonable  temperatures.  Ac¬ 
cordingly,  we  have  developed  and  used  such  an  approach  to  analyze  Hall  and  conduc¬ 
tivity  data  taken  on  samples  grown  by  the  United  States  Army  Electronics  Technology 
and  Devices  Laboratory  (ET&DL),  as  well  as  on  specimens  previously  discussed  by 
Martin. ^ 

EXPERIMENTAL 

The  ET&DL  specimens  were  qrown  by  means  of  a  modified  liquid  (molten  B2O3)  en¬ 
capsulated  Czochralski  technique  that  employs  in  situ  compounding  and  results  in  low 
silicon  content  (typically,  3  X  1014  at  cm-3).^~  The  material  was  cut  into  Greek 
cross  structures  by  ultrasonic  milling,  and  contacted  with  tin, which  was  annealed  at 
450OC  for  5  minutes  in  flowing  hydrogen.  Van  der  Pauw  measurements  of  conductivity 
and  Hall  coefficient  were  made  at  room  temperature  with  a  fully  guarded  system  capa¬ 
ble  of  measuring  conductivities  of  less  than  10-12  Siemens/cm. 6,/ 

THE  BASIC  EQUATIONS 

In  the  absence  of  a  magnetic  field,  the  Hall  coefficient  is  given  by 

R  =  (Rp°p2  +  Rn°n2)/(on  +  °p)2  (1) 

and  the  conductivity  by 

0  =  °n  +  °p  =  newn  +  p%  =  yne  +  pb~^  ^ 

where  e  is  the  electron  charge,  pn  and  pp  are  the  respective  electron  and  hole 
concentrations,  un  and  jpp  their  mobilities,  and  b  their  ratio  un/up.  The  elec¬ 
tron  and  hole  Hall  coefficients,  Rp  and  Rn,  may  be  written  as  -rn/ne  and  rp/pe,  re¬ 
spectively.  The  scattering  factor;  r,  is  very  difficult  to  measure  or  calculate 
precisely.  However,  general  considerations  indicate  that  it  usually  does  not  vary 
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appreciably  from  unity, 3,8-12  n0r  substantially  affect  the  outcome  of  solutions  in 
the  MCR.  For  the  most  general  case,  (rn  *  rp),  the  Hall  concentration  is  given  by 

H  =  rn/Re  =  (n  +  (up/wn)p)2/(p«'1(up/nn)2-n)  =  (n  +  pb'1)2/(pa'i  b'2-n)  (3) 

while  the  Hall  mobility  mh»  defined  as  Ro  is  written  in  the  form 

MH/unrn  =  (pa'1(pp/v‘n)2-n)/(n  +  P^p/yn))  =  (Pa  lfc)  2~n)/(n  +  Pb  ^ 

where  a  =  rn/rp*  Also  interest  is  the  ratio 

o/an  =  1  +  (p/n)(up/un)  =  1  +  b-1p/n  .  (5) 


Away  from  the  band  edges,  the  electron  and  hole  concentrations  vary  with  Fermi  ener¬ 
gy  and  temperature  T  as 


n  =  Nc(exp((Eg  -  Ef)/kT))  1  and 

(6) 

p  =  Nv(exp(Ef/kT))_1 

(7) 

If  the  energy  of  the  valence  band  edge  Ev  is  taken  to  be  zero,  that  of  the  conduc¬ 
tion  band  edge  Ec  is  then  equal  to  the  bandgao  ,  Ec,  which  varies  with  temperature 
according  to  Eg  =  1.519  -  (5.405  X  10*4)  t2/(t  +  204). 13  Nc  and  Nv  are  the  effec¬ 
tive  densities  of  states  at  the  conduction  and  valence  band  edges,  respectively,  and 
are  given  by  Nc  =  8.80  X  1013  T3/2  (I-t/4000)  and  Nv  =  1.66  X  1015  T3/2  cm-3.14  j0 
complete  the  set  of  equations  needed  to  obtain  the  carrier  concentrations  and  mobili¬ 
ties  from  experimental  values  of  conductivity  and  Hall  coefficient,  a  relationship 
between  electron  and  hole  mobility  is  necessary.  At  400K,  a  constant  electron/hole 
mobility  ratio  of 

b  =  yn/Up  =  15  (8) 

is  used  by  Martin, 4  while  at  295K  the  following  empirical  relationship  is  proposed 
by  Look, 3 

b  =  un/up  =  13  +  9X10‘4yn  .  (9) 

It  is  to  be  noted  that  these  ratios  are  not  very  different  for  most  observed  mobil¬ 
ities.  To  facilitate  comparisons  of  the  present  work  with  analyses  made  by  Martin 
and  Look,  we  will  consider  the  case  a  =  rn/rp  =  i  and  use  equations  (8)  and  (9)  for 
analyses  of  measurements  made  at  400K  and  295K,  respectively.  In  addition,  discus¬ 
sion  of  the  basic  parameters  and  equations  employed  will  be  limited  to  values  below 
the  lattice-limited  values  of  electron  mobility,  6000  cm^/Vs  and  8000  cm^/Vs  at  400K 
and  295K,  respectively. 1“3 


CONDUCTIVITY  CRITERIA  FOR  THE  ONSET  OF  MIXED  CONDUCTION  EFFECTS 

A  particularly  illuminating  and  analytically  convenient  basis  for  the  investi¬ 
gation  of  mixed  conduction  effects  is  obtained  by  plotting  the  variables  a,r/|Re|, 
n  and  p  vs.  the  Fermi  energy  Ef  on  a  semi -log  graph,  while  observing  lattice- 
limited  mobility  constraints.  Conductivity  and  r/|Re|  are  plotted  for  several  val¬ 
ues  over  the  range  of  possible  mobilities.  There  is  only  one  r/|Re|  curve  in  Fig¬ 
ures  la  and  lb,  since  the  mobility  ratio  b  is  taken  to  be  a  constant  here.  In 
Figure  2,  where  b  is  a  function  of  nn,  there  is  a  family  of  r/|Re|  curves.  All 
of  the  foregoing  plots  appear  as  straight  lines  over  most  of  their  ranges,  with  any 
departures  from  linearity  (that  is,  departures  of  a  from  on.  Op  and  r/|Re|from 
p,n)  occurring  only  within  the  mixed  conduction  region.  Such  a  graph  affords  a 
simple,  visually  obvious  and  accurate  establishment  of  mixed  conduction  criteria, 
the  boundaries  of  the  MCR  being  the  outer  limits  at  which  significant  deviations 
from  linearity  occur  for  any  curve.  A  deviation  of  1/ ] Re |  from  n  of  about  8%  is 
sufficiently  small  to  serve  as  the  criterion  for  such  a  boundary,  yet  discernible  on 
graphs  of  reasonable  size.  The  deviations  of  o  from  eynn  at  the  boundaries  are 
too  small  to  be  discernible  (see  Figure  lb  for  example).  The  conductivities  asso¬ 
ciated  with  this  criterion  are  10-5  and  2  X  10-9  S/cm  for  temperatures  of  400K  and 
295K,  respectively.  For  conductivities  above  these  values,  equations  2-4  reduce  to 


r*  ^  cr  -  u  en 

n  n 

(2A) 

r/ 1 Re  |  %  n 

(3A) 

,JH  =  %  r,‘n 

(4A) 

and  we  have  a  solution  for  n  in  terms  of  the  measured  quantity  R.  Solutions  for 
Ef  and  p  follow  from  a  consideration  of  Figure  la.  A  glance  at  Figures  la  and  2 
suffices  to  ascertain  that  for  negative  R,  i.e.,  Fermi  energies  above  the  R  =  0 
line,  mixed  conduction  effects  are  small  for  conductivities  above  the  line  NMC  at 
about  6.4X10-5  and  1 . 7X10-9  Siemens/cm  for  temperatures  of  400K  and  295K,  respec¬ 
tively.  As  shown  in  Figures  lb  and  2,  the  NMC  conductivity  is  defined  by  the  inter¬ 
section  of  the  conductivity  curve  for  the  lattice-limited  mobility  and  the  p  =  n 
line.  This  choice  of  criterion  results  in  deviations  of  o  and  from  un  and  un 
of  about  7%  and  of  R  from  r/ne  of  12%,  which  is  sufficiently  accurate  for  most 
analytical  purposes.  It  has  the  additional  advantage  of  providing  easily  recogniz¬ 
able  graphical  features  to  mark  the  onset  of  mixed  conduction.  The  third,  and  per¬ 
haps  most  important,  advantage  afforded  by  this  choice  will  be  discussed  on  page 
10.  The  region  defined  by  the  p  =  n  and  R  =  0  lines  and  labeled  NRMCR,  represents 
the  MCR  for  negative  R.  In  the  next  section  we  will  expand  upon  the  reason  for 
associating  the  departure  from  linearity  with  the  Fermi  level  where  p  =  n.  At  this 
energy,  the  deviation  of  l/|Rej  from  n  is  about  12%.  Table  1  indicates  how  and 
o  differ  from  and  on,  respectively,  at  this  and  other  key  Fermi  energies  for 
various  values  of  b. 

The  conductivity  of  Martin's4  sample  C2  at  400K  was  about  2.4  X  10_5  Siemens/cm. 
Since  this  is  appreciably  above  tne  NMC  line,  mixed  conduction  effects  are  minimal, 
and  the  measured  Hall  mobility  of  3200  cm^/Vs  should  be  a  very  good  approximation  to 
run.  Martin's  estimate  of  ryn  for  this  sample  was  4500-5000  cm2/Vs.  The  long  ar¬ 
rows  on  Figure  la  indicate  the  system  Fermi  energies  (Efc)  determined  by  the  inter¬ 
section  of  the  dotted  (Re)-l  data  line  for  C2(-4.6  X  1010  cm-3)  and  the  r/|Re|  curve 
where  a  scattering  coefficient  r  of  1  has  been  assumed.  The  intersection  of  the 
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Figure  la.  Carrier  Concentrations,  Hall  Concentration  and  Conductivity  vs. 

Fermi,  T=400K. 

Carrier  concentrations  n,p,  Hall  concentration  r/|Re|,  and  conductiv¬ 
ity  o,  as  functions  of  Fermi  level  at  400K.  The  family  of  a  curves  indicates  how 
o  depends  on  electron  mobility  yn.  Measurements  from  Samples  C2  and  Bl,  Table  3, 
are  analyzed. 
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Figure  2.  Carrier  Concentrations,  Hall  Concentration  and  Conductivity  vs. 
Fermi,  T=295K. 

Hall  concentration  and  conductivity  as  functions  of  Fermi  energy  and  the 
parameter  electron  mobility  at  295K.  Measurements  from  sample  Cr2,  Table  2, 
are  analyzed. 
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dotted  conductivity  data  line  (2.4  X  10-5  S/cm)  with  these  system  Fermi  levels  (long 
arrows)  at  points  Q1  and  Q2  determine  the  conductivity  curves,  and  thus  the  electron 
mobilities  un. 

The  conductivity  curve  that  runs  through  Q1  has  an  electron  mobility  of  about 
3200  cm^/Vs,  that  is,  y^,  as  expected.  Interpolation  between  curves  can  be  carried 
out  using  equation  (2)  and  the  carrier  concentrations  implied  by  the  system  Fermi 
levels.  Clearly,  Q2  is  well  above  the  conductivity  curve  defined  by  the  lattice- 
limited  mobility  and  represents  an  inadmissible  solution. 

The  conductivities  of  samples  B1  and  Cr2,  displayed  by  the  heavy  horizontal  dot¬ 
ted  lines  in  Figures  la  and  2,  respectively,  are  below  the  NMC  lines  of  these  fig¬ 
ures  and  will  be  considered  in  the  chapter  on  mixed  conduction  analysis,  page  8. 

For  positive  Hall  coefficients,  that  is,  for  Fermi  energies  less  than  Ef(R=0), 
the  conductivity  criterion  PMC  is  higher  with  values  of  2  X  10"5  and  5  X  10-$ 
Siemens/cm  for  400K  and  295K,  respectively.  The  PMC  line  is  determined  by  the  inter¬ 
section  of  the  conductivity  curve  for  the  lattice-limited  mobility  with  the  "low- 
energy"  border  (LBMCR  of  the  MCR.  The  latter  was  chosen  to  exhibit  the  same  devi¬ 
ation  (^12%)  of  r/ 1 Re  I  from  carrier  concentration  (here  p)  as  obtained  at  the  enerqy 
where  n  =  p  (the  high-energy  border  of  the  MCR).  (See  Figures  lb  and  2.) 

As  may  be  inferred  from  a  review  of  the  data  of  Table  1,  for  positive  '  !  co¬ 

efficients,  ap/a  is  greater  than  90°/  for  b  greater  than  10. 

TABLE  1.  IMPORTANT  CHARACTERISTICS  OF  THE  ANALYSIS 

Relative  magnitudes  which  indicate  the  degree  of  mixed  conduction  at  key  po  •*  ,n 
the  NRMCR  as  a  function  of  un/up.  A  scattering  factor  r  of  1  is  assumed. 
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V1 

16. 6a 

1.06 

0.00 

2.00  0.47 

1  .06 

0.94 

15.0 

1.07 

0.00 

2.00  0.47 

1.06 

0.93 

10. 0b 

1.10 

0.00 

2.00  0.45 

1.10 

0.90 

a.  16.6  corresponds  to  pn  of  4000  in  equation  9. 

b.  Below  minimum  predictable  yR  of  equation  9,  but  interesting  for  perspective. 


Applications  of  Conductivity  Criteria 

The  results  of  room-temperature  measurements  on  ET&DL  samples,  together  with  a 
comparison  of  Look's  old  and  new  mixed  conduction  analyses  with  the  present  analysis, 
are  displayed  in  Table  2.  Samples  U1-U4  are  undoped,  and  both  the  present  criterion 
for  R  0,  (o  >  2  X  10-9  s/cm)  and  Look's  new  analysis  indicate  that  mixed  conduc¬ 
tion  effects  are  minimal,  and  n  and  gn  are  given  by  (eR)*l  and  Ro,  respectively. 
The  results  of  Look's  old  analysis  differ  from  these  conclusions  and  include  predic¬ 
ted  hole  mobilities  well  over  the  lattice-limited  value  of  400  cm3/Vs. 

Look's  new  analysis,3  which  limited  itself  to  room-temperature  measurements  of 
samples  exhibiting  negative  Hall  coefficients,  established  a  minimum  conductivity 
criterion  of  2  to  2.5  X  10*9  S/cm  for  MC  analysis.  This  agrees  quite  well  with  the 
criterion  established  by  the  minimum  visible  deviation  of  r/'Re'  from  n 
(2  to  2.5  X  10*9  S/cm)  in  Figure  2;  the  vertical  dotted  line  defined  by  the  8T  de¬ 
viation  of  the  Hall  curve  from  linearity  results  in  a  conductivity  criterion  of 

2  X  10*9  s/cm.  Look's  technique  employs  a  family  of  curves  relating  to  un 

through  the  parameter  resistivity,  and  is  basid  on  equations  (1),  (2),  and  (9)  and 
a  tentative  value  of  ni  (2.6  X  10*6  cm-3),  determined  as  a  best  fit  to  measurements 
made  at  295K.  The  value  of  n-j  can  be  obtained  either  from  the  intersection  of  lines 
n  and  p  in  Figure  2,  or  from  the  following  expression: 

n.  =  (np ) X/2  =  (NcNv  exp  (-Eg/kT))1^  (10) 

where  the  effective  densities  of  states,  Nc  and  Nv,  of  Blakemore^  are  used  to  ob¬ 
tain  n-j  =  1.25  X  106  cm*3  at  295K.  Thus,  the  resulting  conductivity  criteria  are 
not  strongly  dependent  on  the  exact  choice  of  n-j  nor,  as  can  be  seen  from  Table  1, 
on  the  mobility  ratio  b.  Table  1  shows  that  for  mobility  ratios  b,  as  low  as 

10  pn  is  close  to  uH  and  -n  to  a  at  the  intrinsic  Fermi  level. 


Dependence  of  Carrier  Dominance  on  Fermi  Energy 

The  three  long  solid  lines  in  Figure  lb  mark  energies  defined  in  Table  1.  Be¬ 
cause  the  R  =  0  and  p  =  n  lines  are  well-defined,  easily  recognizable,  and  are  suf¬ 
ficiently  close  to  the  energies  where  discernible  deviations  of  o  from  linearity 
first  occur,  they  are  convenient  borders  for  the  region  where  neither  carrier  domi¬ 
nates  conduction.  These  lines  also  separate  regions  of  decreasing  degrees  of  p-dom- 
inance  arising  with  increasing  Fermi  energy.  In  the  region  where  R  >  0,  complete 
p-dominance  prevails;  that  is,  the  Hall  effect  conductivity  and  carrier  concentra¬ 
tion  are  all  hole-dominated,  and  the  sample  is  denoted  as  being  p-type  of  the  third 
degree,  or  ppp.  In  the  intermediate  region  between  R  =  0  and  the  locus  of  the  con¬ 
ductivity  minima,  R  <  0  but  ap  >  an  and  p  >  n,  so  we  have  p-dominance  of  the  second 
degree,  or  ppri.  The  following  region  has  as  its  upper  bound  the  Fermi  energy  at 
which  the  n  and  p  lines  cross.  Here  we  have  R  <  0  and  op  <  on,  but  p  is  still 
greater  than  n,  so  that  only  first  degree  or  p  concentration  dominance  results 
and  is  denoted  by  pnn.  In  the  region  labeled  nnn,  R  <  0  op  <  an  and  p  <  n;  thus, 
complete  n-dominance  prevails. 


ANALYSIS  OF  MIXED  CONDUCTION  EFFECTS 

For  conductivities  below  the  established  criteria  we  will  graphically  solve 
equations  (2),  (3),  (6),  (7),  and  (8)  or  (9),  simultaneously,  to  resolve  the  elec¬ 
tron  and  hole  concentrations  and  mobilities.  The  whole  process  may  be  considered 
equivalent  to  solving  the  following  quartic  equation  for  n  and  equations  (6),  (7) 
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(11) 


and  (2),  respectively,  for  Ep  p  and  yn: 

(Reb^)n^  -  (rb^)n^  +  (Rebn^)n^  +  (n^r)  n  +  (Ren.^)  =  0 

As  we  shall  see,  some  "programmed  solutions"  to  equation  (11)  can  be  misleading,  and 
it  is  always  best  to  have  a  complete  picture  of  the  processes  involved,  as  afforded 
by  Figures  1  and  2. 


Constant  Mobility  Ratio  (Measurements  Made  at  400K) 

Let  us  consider  the  simplest  case  first,  that  of  equation  (8)  used  at  400K  by 
Martin. ^  According  to  Descartes'  rule,  equation  (11)  has  two  changes  of  sign,  and 
therefore  two  real  solutions.  As  indicated  in  equation  (1 1 ) ,  if  b  is  constant,  n 
and  therefore  Ef,  is  independent  of  un  or  a.  Thus,  there  exists  only  one  curve 
for  R,  as  shown  in  Figure  la,  and  a  measurement  of  R  immediately  yields  the 
double-valued  solution  for  n,  p  (uy/pn»  °/°n  ^rom  equations  (4)  and  (5))  and  the 
system  Fermi  levels.  The  conductivity  curves,  or  equation  (2),  then  yield  the 
double  solution  for  un.  Hall  and  conductivity  data  from  Sample  B1 ,  Table  3,  are 
analyzed  in  Figure  la. 

As  can  be  seen  from  Figure  lb,  the  entire  Hall  curve  to  the  left  of  its  minimum 
lies  in  the  NRMCR  for  negative  values  of  R.  Since  one  of  the  values  for  any  ambig¬ 
uous  solution  for  a  Hall  measurement  must  lie  on  this  portion  of  the  curve,  it  must 
fall  within  the  NRMCR.  Thus,  the  maximum  conductivity  that  can  be  related  through 
a  Fermi  level  to  that  value  of  R  is  defined  by  the  intersection  of  the  R=  0  line 
with  the  conductivity  curve  corresponding  to  the  lattice-limited  mobility.  Since 
line  NMC  denoting  this  conductivity  also  intersects  the  maximum  conductivity  curve 
at  the  p  =  n  line,  the  latter  becomes  the  high  energy  border  of  the  double  solution 
region.  Thus,  NRMCR,  defined  earlier,  also  corresponds  to  the  region  for  which 
there  are  two  real  solutions  to  the  quartic.  Hence,  the  rationale  for  the  previous¬ 
ly  discussed  p  =  n  (or  12%  deviation)  criteria  for  mixed  conduction  becomes  clear. 

Examples  of  such  double-valued  solutions  are  represented  by  the  short  and  long 
arrow  pairs  in  Figure  la  for  Samples  B1  and  C2,  respectively.  In  each  case,  the 
sample  Fermi  level  Ep$  is  found  on  the  r/|Re|  curve  from  the  measured  value  of  the 
Hall  coefficient.  The  conductivity  curve  passing  through  the  conjunction  of  Eps 
with  the  measured  conductivity  then  corresponds  to  the  yn  of  the  sample.  Points  P 
indicate  the  double-valued  solution  of  yn  for  Sample  Bl,  namely,  yn  =  2100  and 
5200  cm^/Vs,  which  satisfy  the  measured  conductivity  2.7  X  10“6  at  tne  two  Fermi 
energies  implied  by  the  measurement  of  1/ ( Re j .  Similarly,  points  Q  indicate  the 
solutions  for  specimen  C2,  but,  in  this  case,  the  solution  corresponding  to  point 
Q2  would  be  well  above  the  lattice-limited  mobility  and  is,  therefore,  not  admis¬ 
sible.  Consideration  of  point  Q2  is,  of  course,  also  ruled  out  by  the  conductivity 
criterion  established  above. 

For  positive  Hall  coefficients,  two  admissible  solutions  occur  in  the  PRMCR 
for  conductivities  below  line  1.1X10”^  S/cm  in  Figurel(where  PDV  is  an  abbreviation 
for  positive  double-value).  However,  significant  MC  effects  on  the  approximation  of 
r/ | Re |  to  p  continue  to  exist  for  conductivities  as  large  as  2  X  10“5  S/cm  or  out 

to  the  Fermi  energy  labeled  LBMCR.  Therefore,  the  PDV  line  does  not  intersect  the 

conductivity  curve  corresponding  to  the  lattice-limited  mobility  at  LBMCR,  the  lower 

border  of  both  PRMCR  and  MCR.  This  means  that  the  criterion  for  the  onset  of  mixed 

conduction  effects  does  not  coincide  with  the  criterion  for  the  onset  of  double¬ 
valued  solutions  as  it  does  for  negative  Hall  coefficients. 
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Three  analyses  of  Martin's  measurements  will  be  discussed,  and  the  results, 
principally  carrier  concentration,  will  be  listed  in  Tables  3  through  5.  All  involve 
the  assumption  that  un/uB  =  15  and  r  =  1.  The  first  analysis  is  based  on  Figure  1, 
the  second  on  solution  of  equation  (11)  for  the  n-j  (6.22  X  10^  cm_3)  generated  in 
Figure  1,  and  the  third  on  solution  of  equations  (2),  (3),  and  (8)  and  use  of  Mar¬ 
tin's  estimate  of  ryn,  which  he  refers  to  as  "electron  mobility."  The  first  two 
analyses  should,  but  do  not  always,  result  in  fair  agreement.  Unlike  the  first  two, 
the  third  analysis  results  in  only  one  solution  for  each  set  of  measurements  and  a 
large  spread  in  the  n-j  values.  The  average  n^  is,  however,  in  good  agreement  with 
the  values  deduced  or  assumed  in  the  first  two  analyses. 

Table  3  lists  the  two  sets  of  solutions,  including  Fermi  level  and  carrier  con¬ 
centrations  implied  by  each  individual  Hall  measurement.  Figure  1,  and  the  scatter¬ 
ing  coefficient  r.  The  mobilities  are  then  found  from  the  conductivity  measurement 
through  use  of  interpolation  in  Figure  1  or  equation  (2)  and  the  carrier  concentra¬ 
tions. 

The  upper  and  lower  sections  of  Table  1  reflect  scattering  factors  of  1  and 
1.26,  respectively.  The  former  is  the  convenient  approximation  typicallyl_3  made, 
but  is  inconsistent  with  the  measurements  of  samples  A9,  B2,  and  the  r/|Re|  minimum 
of  Figure  1.  The  latter  both  satisfies  these  data  and  is  near  the  center  of  the  es- 
timatedS  range  in  r,  1.2  ±  20".  The  solutions  obtained  for  samples  B1  and  B3  for  an 
r  of  1.26  are  listed  to  indicate  the  size  of  the  changes  in  concentrations  and  mo¬ 
bilities  resulting  from  this  increase  in  scattering  factor.  The  change  for  solu¬ 
tions  inside  the  mixed  conduction  region  is  small,  and,  of  course,  this  is  the  pri¬ 
mary  area  of  interest  here.  Outside  the  mixed  conduction  region,  uncertainty  in  the 
scattering  factor  has  a  much  greater  effect  on  the  estimation  of  un  from  the  ryn  or 
UH  data,  such  as  Martin  obtains  from  epilayer  measurements. 

Table  4  lists  calculator-generated  solutions  of  equation  (11)  based  on  the 
method  of  steepest  descent. ^  Four  solutions  are  generated,  all  are  complex  for 
1/ | Re |  <  4.0  X  10^;  for  1/ | Re |  >  4.1  X  10^  cm_3,  two  are  real  and  two  are  complex 
with  comparable  real  and  imaginary  components.  The  amplitudes  of  the  complex  solu¬ 
tions  are  the  same,  unlike  those  of  the  real  solutions.  As  the  value  of  l/;Re|  ap¬ 
proaches  and  drops  below  4.0  X  10^  cm'3,  the  two  real  solutions  merge  and  two  com¬ 
plex  solutions  having  the  same  real  part  emerge.  As  Figure  1  correctly  indicates, 
this  merger  should  occur  as  l/|Pej  drops  below  the  minimum  value  of  5.06.  This  de¬ 
lay  in  "coupling”  the  real  roots  generated  by  the  computer  results  in  two  real  solu¬ 
tions  for  Sample  B2  which  fail  to  satisfy  equation  3  for  R. 

Table  5  lists  Martin's  measurements  and  estimates  of  ryn,  as  well  as  the  re¬ 
sulting  carrier  concentrations  and  range  of  n-j  determined  for  r  =  1  from  equations 
(2),  (3),  and  (8).  Using  a  knowledge  of  the  concentrations  of  dominant  dopants  in 
his  bulk  material,  together  with  a  plot  of  electron-free  carrier  concentration  vs. 
Hall  mobility  in  very  high-purity  VPE  samples,  Martin  estimated  the  electron  mobili¬ 
ty  in  his  bulk  samples.  Samples  A9  and  A10  are  dominated  by  deep  donor-level  (re¬ 
ferred  to  as  EL2)  concentrations  of  a  few  1016  cm-3,  and  were  estimated  to  have  elec¬ 
tron  mobilities  of  5000  cm^/V.  Samples  B1  through  B5  were  dominated  by  chromium 
concentrations  near  10^  cm_3,  and  were  estimated  to  have  electron  mobilities  of 
4000  cm2/Vs.  These  estimates  permitted  deduction  of  the  ratio  un/ryn  which,  in 
turn,  was  used  with  b  to  deduce  values  for  the  concentration  ratios  [EL2]/  (N/\-N[j) 
and  [Cr]/  ([EL2]  +  Na-Nd). 

In  a  later  paper, ^  Martin  referenced  results  of  Look's  old  mixed  conduction 
analysis  as  evidence  that  ryn  varies  only  slightly  as  a  function  of  ionized  im¬ 
purity  scattering,  and  as  a  rationale  for  treating  ryn  as  constant  at  about  4000 
cm^/Vs.  As  shown  in  Table  5,  the  estimations  from  reference  4  imply  a  variable  n-j 
ranging  from  2.9  to  11  X  10^  cm"3.  Interestingly,  the  average  n-j  is  6.1,  which 


TABLE  3.  ANALYSIS  OF  MARTIN'S  MEASUREMENTS  USING  FIGURE  1 


Carrier  concent  rat  ions  and  relevant  ratios  Implied  by  Fig.  1  and  Hall 
effect  measurements  from  Ref.  4;  and  mobility  determinations  implied  by  these  carrier 
concentrations  and  the  associated  conductivity  measurements.  Scattering  factors,  rt 
of  1.0  and  1.26  are  used. 


c.  LEC-grown  crystal,  heavily  doped  with  chromium 

d.  Of  greater  interest  here  are  the  U„/w  and  of  o  ratios  of  0.88  and  1.01, 

H  p  p 

respectively 

e.  Absolute  value  less  than  5.06  cm  which  is  the  minimum  allowable  value  of 
1 / | Re |  for  parameters  of  Fig.  1 

f.  Corresponds  to  minimum  value  of  1.264/ | Re |  for  parameters  of  Fig.  1 

g.  Larger  than  lattice-limited  mobilities 
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TABLE  4.  ANALYSIS  OF  MARTIN’S  MEASUREMENTS  USING  COMPUTER-GENERATED 

SOLUTION  TO  QUARTIC  EQUATION  11 

Q 

Calculator-generated  solutions  to  equation  11  for  r  =  1  and  n^  =  6.244  X  10  cm  . 
(Compare  with  electron  concentrations  in  Table  3. 

SAMPLE  (eR)'1  Real  n 

(10^cm~3)  (lO^cm-^) 


A9 


-4 


a.  Complex 

b.  Complex 


A10 


-11 


a.  10.7 

b.  1.05 


B1 


-8.5 


a.  8.15 

b.  1.18 


B2 


-4.5 


a.  3.57 

b.  1.87 


B3  -17  a.  16.8 

b.  0.89 

B4  -28  a.  27.9 

b.  0.76 

B5  +312  a.  0.14’ 

b.  0.32 


i.  The  corresponding  value  of  p  is  278  x  lO^cm"^. 
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TABLE  5.  ANALYSIS  OF  MARTIN'S  MEASUREMENTS  USING  MARTIN'S  ESTIMATES 

OF  ELECTRON  MOBILITY 
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Figure  3.  Carrier  Concentrations,  Conductivity  and  Conductivity 
Product  vs.  Fermi  Energy  as  a  Function  of  Electron  Mobility, 
T=295K. 

Measurement  product,  aun/r  or  o2R/r  and  conductivity  as  functions 
of  Fermi  energy,  and  the  parameter  electron  mobility  at  295K. 
Measurements  from  sample  Cr2  are  analyzed. 


agrees  well  with  the  value  6.2  deduced  from  Figure  1.  For  three  cases  in  Table  3, 
B3  and  B5  (for  r  =  1)  and  B2  (for  r  =  1.264),  the  predicted  values  of  run  (which 
are  in  the  MCR)  agree  reasonably  well  with  the  estimates  made  by  Martin, 4  and  fluc¬ 
tuations  in  such  estimates  are  to  be  expected.  Clearly,  the  establishment  of  a  val¬ 
ue  for  a  fundamental  quantity  like  the  intrinsic  carrier  concentration  at  the  tem¬ 
perature  of  interest  is  quite  important,  and  the  present  analysis  affords  a  useful 
vehicle  for  such  considerations.  (See,  also,  Blakemore.^)  If  we  choose  to  assume 
values  for  the  intrinsic  carrier  concentration  and  b,  then  the  measurements  of 
R  and  a  in  Table  3  can  be  satisfied  by  two  plausible  solutions,  both  of  which 
must  be  considered.  In  each  case,  one  solution  is  outside  the  mixed  conduction  re¬ 
gion,  except  for  analyses  of  B2  and  A 9  for  r  =  1.26.  In  one  case,  A10,  only  the 
outer  solution  has  a  mobility  below  the  lattice-limited  value.  We  will  consider 
later  other  techniques  to  permit  discrimination  between  the  ambiguous  solutions. 


Variable  Mobility  Ratio  b(un)  (Measurements  Made  at  295K) 

The  above  analysis  of  Figure  1  assumed  a  constant  mobility  ratio,  b  =  15,  and 
produced  a  single  curve  based  on  the  Hall  coefficient  which,  in  turn,  depended  only 
on  Ef.  In  contrast,  this  analysis  assumes  a  variable  b(pn)  of  equation  (9)  to 
produce  a  family  of  such  Hall  curves  which  depend  on  yn  or  a  as  well  as  on  £f. 
Figure  2  displays  three  members  of  this  family  of  Hall  curves  for  electron  mobili¬ 
ties  of  2000,  4000  and  8000  cm^/Vs.  The  range  in  Fermi  energy  for  this  plot  had  to 
be  reduced  to  half  that  of  the  other  figures  to  permit  clear  display  of  three  mem¬ 
bers  of  the  family.  For  this  reason,  the  c2|R|  plots  of  Figure  3  may  be  preferred 
for  analysis.  Unlike  that  for  1 / ( Re [ ,  these  plots  have  a  symmetry  similar  to  that 
of  the  conductivity  curves,  and  embody  what  is  perhaps  the  simplest  form  for  the  re¬ 
lationship  between  the  measured  Hall  coefficient  and  the  transport  properties 

Ra2/r  =  yHo/r  =  eun2  (pb"2-n)  .  (12) 

Either  Figure  2  or  Figure  3  can  be  used  for  mixed  conduction  analysis  for  the  mobil¬ 
ity  ratio  b  of  equation  (9). 

Figures  2  and  3  depict  analyses  of  the  conductivity  and  Hall  measurements  of 
sample  Cr  2,  Table  2.  The  solution  entails  determination  of  a  mobility  which  allows 
satisfaction  of  both  the  conductivity  and  Hall  coefficient  measurements  at  the  same 
Fermi  level.  The  analysis  in  Figure  2  is  simplified  by  the  existence  of  a  single 
r/ | Re |  curve  for  energies  above  the  MCR.  Thus,  the  Fermi  energy  of  the  solution 
outside  the  MCR  (single  arrow)  is  uniquely  determined  by  the  Hall  measurement;  the 
mobility  (about  2000  cm^/Vs)  is  then  implied  by  the  conductivity.  The  broad  arrow 
in  the  MCR  indicates  the  ambiguity  in  the  Fermi  energy  as  determined  solely  from  the 
r/ 1 Re  1  curves  and  measurement.  It  is  immediately  clear,  however,  from  this  broad 
arrow  and  the  associated  conductivity  measurement,  that  the  solution  in  the  MCR  im¬ 
plies  mobilities  above  the  lattice-limited  values.  Thus,  even  when  a  conductivity 
measurement  implies  an  ambiguous  solution,  it  is  often  rendered  unique  by  this  mo¬ 
bility  consideration. 

The  technique  for  resolving  an  ambiguous,  admissible  solution  (Ef  and  yn)  in 
the  MCR  of  Figure  2  is  similar  to  that  shown  in  Figure  3  for  identifying  a  solution 
in  any  region.  The  dotted  lines  in  Figure  3  illustrate  this  process  for  the  same 
data  used  in  Figure  2.  These  lines  connect  each  a^R  curve  at  the  measured  value 
with  thecurve  of  corresponding  mobility  at  the  measured  conductivity.  The  vertical 
line  A-B  corresponds  to  the  Fermi  level  associated  with  the  solution.  As  with 
Figure  2,  this  analysis  results  in  one  solution  which  satisfies  the  measurement  and 
mobility  limitations.  When  there  are  two  admissible  solutions,  one  can  turn  to 
charge  balance  plots,  thermopower,  or  other  measurements  described  by  Look^  to 
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discriminate  between  solutions. 


Look's  new  analysis  admits  of  two  solutions  for  each  of  the  Cr-doped  samples, 

Cr  1  and  Cr  2,  Table  2.  In  each  of  these  cases,  Figures  2  and  3  admit  of  only  one 
solution  and  it  is  outside  the  mixed  conduction  region.  The  inconsistency  in  the 
results  can  only  be  attributed  to  the  difference  in  nis  chosen  by  Look  as  a  best  fit 
to  his  data  and  the  value  implied  by  Blakemore's^  Ny.  Nc  values. 


CHARGE  BALANCE  ANALYSIS 

The  Fermi  level  of  any  semiconductor  is  fixed  by  the  necessity  for  charge  con¬ 
servation;  that  is,  the  number  of  electrons  released  from  donors  and  the  valence 
band  must  equal  the  number  received  by  acceptors  and  the  conduction  band.  This  con 
dition  is  expressed  by: 

lN‘+n=£Nj+p  (a) 


where  N/\  and  Np  are  the  numbers  of  ionized  acceptors  and  donors,  respectively, 
which  are  given  by: 


+  ( NA’  V 

<YV  '  T*T9079\T^xpTnrAT^^l^T7kT) 


(b) 


Here  we  rely  on  Martin's  work  for  the  characterization  of  the  Cr  and  EL2  levels.  In 
particular,  energy  levels  of  0.730  eV  and  0.705  eV  (measured  from  the  valence  band) 
are  adopted  for  Cr  and  EL2,  respectively.  Martin  found  these  values  to  be  consis¬ 
tent  with  g0  gi  values  of  0.93  and  1,  respectively. 

Condition  (a)  can  be  expressed  graphically  by  plotting  both  sides  of  (a)  against 
Ef  on  a  semi -log  plot,  as  done  for  Sample  B3  in  Figure  4.  The  Fermi  level  is  deter¬ 
mined  by  the  intersection  of  (sN/(  +  n)  and  (eNq  +  p).  This  technique  was  introduced 
by  Shockleyl?  and  has  been  used  recently  by  others. 18-20  The  conductivity  curve  for 
the  deep  Fermi  level  solution  of  Sample  B3  and  the  Hall  concentration  curve  are  also 
drawn  in  on  Figure  4  to  permit  comparison  of  the  Fermi  levels  obtained  using  Figure 
la  with  that  using  Shockley  diagram  analysis. 

Such  comparisons  permit  resolution  of  ambiguous  solutions  even  when  impurity 
concentrations  are  known  only  approximately.  For  example,  specimen  B3  was  lightly 
doped  with  Si  to  a  concentration  of  about  3  X  1016  cm--*  and  was  estimated  to  have  a 
similar  EL2  concentration;  B3  was  doped  to  a  Cr  concentration  of  3  X  10*7  cm-3. 

Since  this  is  about  30  times  greater  than  the  total  of  all  other  acceptors,  the  point 
of  intersection  of  the  +  n)  curve  and  the  (iNj  +  p)  curve  is  little  affected  by 
the  inaccuracy  of  our  estimate  for  the  total  shallow  acceptor  concentration  for  Fermi 
energies  of  interest  here,  e.g.,  those  above  the  discontinuity  in  the  Hall  curve. 

Thus  we  obtain  the  estimate  shown  of  the  Fermi  level  as  determined  from  the  impurity 
concentration  information.  The  deep  and  shallow  Fermi  levels  labeled  Ef (H)  and  E'f 
are  determined  by  the  intersections  A  and  A',  respectively,  of  the  Hall  measure¬ 
ment  with  the  Hall  curve.  The  relatively  small  displacement  (0.005  eV)  of  the  Fermi 
level,  determined  from  impurity  concentrations  from  the  deep  level  solution  at  A, 
resolves  the  ambiguity  in  the  dual  value  solution. 


Figure  4.  Charge  Balance  Plot 


Charge  balance  plot  representing  the  ionized  impurity  structure 
of  Sample  B3.  Also  shown  are  the  conductivity  curve  representing  the  mo¬ 
bility  associated  with  the  deep  Fermi  level  solution  found  using  Figure  la 
and  the  Hall  concentration  curve. 
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Conversely,  if  the  Fermi  level  is  known  from  the  conductivity  and  Hall  concen¬ 
tration  analysis  of  Figure  1  (using  other3  electrical  measurements  to  resolve  any 
ambiguities),  and  if  all  the  donor  and  the  shallow  acceptor  concentrations  are  known, 
then  the  Cr  concentration  could  be  deduced  using  Shockley  diagrams.  More  generally, 
the  ramifications  of  changes  in  the  relative  impurity  concentrations  are  subject  to 
immediate  visual  inspection. 


SUMMARY 

We  present  an  analysis  of  mixed  conduction  effects  based  on  energy  gap,  recent¬ 
ly  deduced  electron-to-hole  mobility  ratios  used  in  related  analyses,  and  newly  re¬ 
vised  calculations  on  the  effective  densities  of  states.  Results  of  each  of  the 
analyses  on  appropriate  data  sets  are  compared.  The  technique  presented  permits 
consideration  of  charge-balance  effects  and  impurity  concentrations  on  the  same 
graph,  as  well  as  the  establishment  of  simple  criteria  for  the  necessity  of  mixed 
conduction  analysis.  This  affords  a  means  for  immediate  visual  discrimination  be¬ 
tween  double-valued  solutions  of  carrier  concentrations,  or  visual  consideration  of 
the  impurity  structure  dominating  compensation  in  the  material. 
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